Collective modes of an anisotropic quark-gluon plasma induced by relativistic jets 
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We discuss the characteristics of collective modes induced by relativistic jets in an anisotropic 
quark-gluon plasma(AQGP). Assuming a tsunami- like initial jet distribution, it is found that the 
dispersion relations for both the stable and unstable modes are modified substantially due to the 
passage of jet compared to the case when there is no jet. It has also been shown that the growth 
rate of instability first increases compared to the no jet case and then completely turned into 
damping except the case when the jet velocity is perpendicular to the wave vector in which case the 
instability always grows. Thus, the introduction of the jet in the AQGP, in general, might to faster 
isotropization for the special case when the wave vector is parallel to the anisotropy axis. 
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I. INTRODUCTION 

The primary goal of the ultra-relativistic heavy-ion col- 
lision experiments at BNL RHIC and at CERN LHC is 
to study the properties of a deconfined state of the QCD 
matter, commonly known as quark-gluon plasma(QGP). 
According to lattice calculation, the novel state of matter 
is expected to be formed when the temperature of the nu- 
clear matter is raised to the critical value T c ~ 170MeV, 
or the energy density of the nuclear matter is raised to 
above 1 GeV//m 3 . High energy partons behave as hard 
probes which are produced in the early stage of the colli- 
sion due to hard scattering. In relativistic heavy-ion colli- 
sion jets with high transverse momentum travels through 
the hot and dense medium and it losses energy by col- 
lisional( interaction with thermal quark and gluon) and 
radiative processes (bremsstrahlung). This phenomena is 
commonly known as jet quenching, because in the direc- 
tion of propagation of the jet one observes a decrease of 
high energy hadrons and increase in the number of soft 
hadrons. In addition, the passage of jet influences the 
collective modes of the system as we shall sec below. 

In Refs (H-dli the effect of a jet of particle on the 
plasma properties, propagated through the equilibrated 
and isotropic QGP has been analyzed, where initial con- 
figuration of jet is assumed to be colorless, electrically 
neutral and described by tsunami-like distribution Q in 
momentum space. Non-equilibrium jet of particles while 
traveling through QGP destabilizes the plasma produc- 
ing the collective gauge fields. These gauge modes might 
have several branches and some of these might even be 
unstable. The most important among those are the 
modes which grow exponentially in time. When the in- 
stability occurs, the kinetic energy of the particles is con- 
verted to the field energy which speeds up the equilibra- 
tion process and also leads to faster isotropization [f|-|8| of 
the QGP. The main assumption for such phenomenon to 
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occur is that the time to generate the growth of the gauge 
fields is smaller than the hadronization time. Plasma in- 
stabilities fully develop on the time scale of the order of 
t ~ (6.7-12.5)/w t 0,||, where uj t is the total angular fre- 
quency of the whole system. In a weak coupling regime, 
the time scale is 1 - 2 fm/c at T - 350MeV. 

In the present work we shall concentrate on the collec- 
tive modes in AQGP induced by relativistic jets. In the 
early stage of heavy-ion collision due to rapid longitudi- 
nal expansion, the plasma cools faster in the longitudi- 
nal direction leading to (p 2 L ) « (pt)- Such momentum- 
space anisotropy leads to collective modes having charac- 
teristic behavior distinct from what happens in isotropic 
plasma which has been extensively studied in [j| 
where it is shown that the gluonic collective modes can be 
unstable. When a relativistic stream of particles interact 
with the non-equilibrium plasma, with an anisotropic dis- 
tribution in momentum space, the behavior of collective 
modes change as will be demonstrated in the following. 

In studying the evolution of such a system, we use the 
method of the plasma physics within the frame work of 
the quark-gluon transport theory 0, in weak coupling 
regime, i.e. g << 1. The time scale for the evolution of 
collective modes is assumed to be much shorter than the 
inter-particle collision time. In this approach we have 
neglected hard mode interactions, assuming that the in- 
teractions between jet and plasma is only mediated by 
mean gauge fields. 

Kinetic instability can occur due to the interaction 
of plasma and jet partons. They are initiated either 
by charge or current fluctuation. In the first case, the 
electric field is longitudinal i.e. the field is parallel to 
the wave vector k (E||k), while in the second case the 
field is perpendicular to k (E _L k). For this rea- 
son the corresponding instabilities are called longitudi- 
nal and transverse instabilities respectively. Since the 
electric field plays a crucial role in the generation of 
longitudinal modes, they are also called electric modes, 
while the transverse modes are called magnetic modes. 
The magnetic mode known as filamentation- or Weibel- 
instability (HJ appears to be relevant for the QGP 0, H| . 
In momentum-space anisotropic plasma, the growth rate 
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of magnetic instability is maximum in the direction of 
anisotropy • So we concentrate on this special case 
in which the momentum of the collective mode is in the 
direction of the anisotropy. We expect that the disper- 
sion relations will be affected due to the passage of jet 
in AQGP. In fact, we shall see that the growth rate of 
unstable modes first increases compared to the case of no 
jet scenario and then becomes damped. 

The organization of the paper is as follows. In section 
2 we briefly recapitulate Vlasov-type transport equation 
for the colored particles. In section 3 and 4 we shall show 
how the dispersion relations are modified in the presence 
of a tsunami like jet. Section 5 will be devoted to discuss 
the results followed by summary in section 6. 



II. VLASOV-TYPE TRANSPORT EQUATION 
FOR THE COLORED PARTICLES 



Transport theory provides a natural framework to 
study equilibrium and non equilibrium plasmas. We 
review the transport equations obeyed by the distribu- 
tion of colored particles and the classical fields. The 
(anti)quark distribution functions f{p,x) and f(p,x), 
which are iV c x N c hermitian matrices, belong to the fun- 
damental representation of the SU(N C ) group. The glu- 
ons are the adjoint representation of SU(N C ), and their 
distribution function Q{p,x) is a — 1) x (N% — 1) 
matrix in color space. The partials of jet consisting 
only quarks, and the corresponding distribution function 
can be described by a N c x N c matrix in color space as 

Wjetfax). 

In the fundamental representation, the color four- 
current j M is expresses as 

3^{X) = ~ ~ I Af{P:X)-f{jp, X ) 



N, 



■Tr[f{p,x) ~ f(p, %)} 



2r a Tr[T a g{p,x)} 



(1) 



where r°, T a with a = 1, ...N 2 - 1 are the SU(N C ) group 
generator in the fundamental and adjoint representations 
with Tr[r a r b ] = \8 ab , Tr[T a T b ] = N c 5 ab and g is the 
QCD coupling constant. Here we use the notation 

d A p 



(2tt> 



:2Q(p )8(p z ) 



(2) 



The transport equations of the quarks, antiquarks, glu- 
ons and the jet particle read as [ll| 

p»D^(p,x) + 9 -p»{F^(x),d;f{p,x)} = C, (3) 

p^Djip, x) ~ ^{F^ix), d%f(p, x)} = C, (4) 

p^v^g(p,x) + 9 -p»{F^{x),d;g{ P ,x)} = c g , (5) 

p^D^W je t(p, x) + ^{F^(x),d;W jet (p, x) = CV (6) 



where T^v = d^A v — d^A^ — ig[A^,, A u ] is the gluon field 
strength tensor, while F^ v represents the strength ten- 
sor in the fundamental representation and the covariant 



derivatives D M and 2? M act as 



ig[Afj,{x), ...], 

with gauge field = A^{x)r a and A^ = A£T a . 

Now, we will study the Vlasov approximation, or colli- 
sionless dynamics of the color conductor where the mean 
free time of the collisional processes are lager then the 
time scale of evolution of the whole system. Therefore 
the collision terms C, C, C a and Cw are neglected. The 
characteristic time scale [13[ of the system's evolution 
due to interparton collision is t} lar d ~ (g 4 T 'ln(l / g)) _1 
and t S oft ~ (<7 2 T ln(l / 'g))^ 1 . So our current approxima- 
tions are really only valid for time intervals shorter then 

tsoft- 

In this work we will also assume that initially the 
plasma is colorless and in thermally equilibrated state. 
Considering the distribution functions are deviation from 
their equilibrium states we have, 



f(p, x) = n fd (p) + Sf(p, x), 
I(P, x) = %d(p) + 5f(p, x), 
G(p,x) = n be (p) + 5G(p,x), 



where 



(P) 



1 



(7) 



(8) 



eP/ T ± 1 

are the Fermi-Dirac and Bosc-Einstcin equilibrium dis- 
tribution functions. 

We also consider color fluctuation of the initial color- 
less distribution function of the jet. Thus 



Wj e t = /jet(p) + SW jet (p,x) 



(9) 



where fjetip) is the initial jet distribution function which 
we will consider to be of colorless tsunami-like form H 



f jet (p)=nu°S^(p-An). 



(10) 



Here n is a parameter proportional to the density; and 
w' 1 is the four-velocity. The parameter A fixes the scale 
of energy of particles. 

With the linear approximation of transport equa- 
tions (Vlasov approximation) one can solve the polariza- 
tion tensor for particles species a 0, H, 0, HH : 



n a (k)-g I p [g p k + i 



(11) 



where a specify the quarks, antiquarks, gluons or partials 
of jet. This tensor is symmetric, H^ u {k) = H vfJ, (k), and 
transverse, k^U^ v = 0. 



III. COLLECTIVE MODES OF ANISOTROPIC 
QUARK-GLUON PLASMA AND THE JET 

To include the local anisotropy in the plasma, one has 
to calculate the gluon polarization tensor incorporating 
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anisotropic distribution function of the particles. This 
subsequently can be used to construct HTL corrected 
gluon propagator which, in general, assumes very com- 
plicated form. Such an HTL propagator was first derived 
in Q in the temporal-axial gauge. 

The spacelike component of the self-energy tensor can 
be written as 



WUk) 



d 3 p 



v j k l 



(12) 



The phase-space distribution is assumed to be given by 
the following ansatz [t| [l(| : 



/(p) = /«(p) = iv(o/i,o(Vp 2 +e(p-fi) 2 )- 



(13) 



Here fi SO is an arbitrary isotropic distribution function. 
N(£) is the normalization constant which is equal to 
VI + £, n is the direction of anisotropy. The parameter 
£ is the degree of anisotropy parameter (— 1 < £ < oo) 

and is given by £ = — 1- Making a change of vari- 

able (p 2 = p 2 (l + £(p.n) 2 )) the spatial components can 
be written as 



where 



dVl V 1 + £(v.n)n' 
(4tt) (l+£(v.n) 2 ) 2 



>dfi SO (p 2 



v j k l 



(k.v + it 
(14) 



= "2^/ d PP 



dp 



(15) 



is the isotropic Debye mass which depends on fi SO , 
The self energy, apart from four-momentum^^ 1 ), also de- 
pends on the anisotropic vector (tt/ 1 = (l,n)) and 
can be cast in a suitable tensorial basis appropriate for 
anisotropic plasma in a co- variant gauge in the following 
way [UGl]: 



where 



m (k) = aA ij + /3B 13 + -/C ij + SD*- 



A ij = gU /,-'/,-'/k 2 , 

B ij = FF/k 2 , 

C ij = h l ti 3 /fi 2 , 

D ij = k l h 3 +k j h\ 



(16) 



(17) 
(18) 
(19) 
(20) 



function 7 and S vanish and a and (3 are directly related 
to the transverse and longitudinal components of the po- 
larization tensor of the plasma respectively. 

The dispersion law for the collective modes of 
anisotropic plasma in temporal axial gauge can be de- 
termined by finding the poles of propagator A y 

A' y '(fc) = „ (22) 

[(k 2 -u 2 )5v - +U l p J (k)} 



Substituting Eqs. (|17H20[) in the above equation and per- 
forming the inverse formula [9j one finds 



A(k) = A A [A-C] + A G [(k 2 



uj + a 



7)B 
-<JD1 



(23) 



The dispersion relation for the gluonic modes in 
anisotropic plasma is given by the zeros of 

A A \k) = k 2 -u 2 + a = 0, (24) 
Ag 1 ^) = (k 2 -uj 2 + a + ~f){/3 -uj 2 ) - k 2 h 2 5 2 = 0. 

(25) 



If we examine the propagators ([24]) and (|25|) in the static 
limit (w — > 0), we find that there are three mass scales [|[: 



0, mi 







m a and m±. In isotropic limit, £ 
and m 2 + = m 2 D . 

The solutions of the above two equations depend on 
mo, u>, k, £ and k.n = cos0 n . For £ > one finds that 
there are at most three stable and two unstable modes 
which depend on 9 n and for £ < 0, there are three stable 
modes but only one is unstable @, [l(| • 

Now, we shall briefly recall the calculation of the po- 
larization tensor induced by the tsunami-like momentum 
distribution of the jet given by Eq. ([T0|) . Substituting 
Eq. (flO|) in Eq.lfTTj) one deduces the following expression 
of the polarization tensor for the jet partons: 



nj et (fc) = -< U 8J + 



k 3 v 



jei 



{uj 2 - k 2 )v 



jet V jet 



k.V 



k.v jet ) 2 



(26) 



where Vj et is the velocity of jet and uj 2 et = ^ s * nc 
plasma frequency of the jet. The dispersion laws of the 
collective modes of the system due to jet are determined 
by searching the poles of the propagator of Eg. (|2"21 by 
replacing IK 3 with H^ et i.e. by finding the solution uj(k). 



where fi l = A 13 n 3 which obeys n.k = 0. 
Now a, /3, 7 and 5 are determined by the following con- 
tractions: 



IV. COLLECTIVE MODES OF THE 
COMPOSITE SYSTEM 



k l lF 3 k J 
//'ll' /. 
/VII' V 
Trrf J 



k 2 /3, 
n 2 k 2 (S 

c;2 



n (a + 7), 
2a + [3 + 7 



(21) 



The expressions for structure functions have been given 
in Ref. @. In the isotropic limit, £ 0, the structure 



Now we study the collective modes of the system due 
to propagation of an energetic jet in an anisotropic quark 
gluon plasma. In very short time regime where the 
Vlasov approximation is valid, the total polarization of 
the system is given by the sum of the two polarization 
tensors: 



nr(fc) = nnfc) 



n^ t (fc) 



(27) 
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The dispersion relation of the collective modes of the to- 
tal system can be determined by solving the equation 

det[(fc 2 - uj 2 )8 13 - fc'V + Hf(fc)] = detlA-^AOP = 0. 

(28) 

The solution of the above equation depends on 
|k|, |vj e t|, k.v jet = cos9 jet , m D , £,k.n = cos6>„ and 
also on 



where 



jet 



(29) 



(30) 



1.5- 



0.5- 

0, 



5 = 10 



1 I 1 I/' / 

■ •/ / . 
.'/ 
■•/ 

.'/ / - 
.•v 
■ •/ 
■'/ 



_l I I I L. 



-i 1 1 1 1 p r 



without jet 
with jet 



_l I I I I I L. 



0.5 1 1.5 2 0.5 1 1.5 2 
k/m k/m 



Now the polarization tensor due to of jet partons can be 
decomposed in the following way: 



where 



n;,;., = n'.v + fir 



(31) 



FIG. 1: (Color online) The dispersion relation for the stable 
A-mode for an anisotropic plasma with jet for different values 
of the anisotropy parameter £ = {1, 10}, r\ = 0.2, 8j et = 7r/2, 
vjet — 0.7 and 8„ — 



a 



"jet 



2(CJ - kVj e t COsOjet) 2 

x [2lo 2 + v 2 et (k 2 ~u 2 ) + Vj et cos 1 

(Vj e t COS 9 jet{k 2 + U) 2 ) — 4kuj)j 
22 V% t COS 2 jet - 1 



jet 



(32) 



J U) — KVjet COS j e t 

Therefore A _1 (fc) in terms of the tensorial basis becomes 

A -1 (fc) = (\L 2 -Lu 2 + a + a')A+(l3 + l3'-uj 2 )B + ~fC + 5'D. 

(33) 

Thus, we obtain an expression for the effective propaga- 
tor as, 

A(fc) = A A A + (k 2 - oj 2 + a + a' + 7)A G B 

+ p + /?'-w 2 )A G - A A }C-5A G T> (34) 



with 



A^ 1 



A 



k 2 -uj 2 + a + a', (35) 
- (k 2 - u? + a + a' + j)(p + 13' - to 2 ) 
- k 2 n 2 S 2 . (36) 



In the following subsection, we will analyze the nature 
of the modes of the whole system. 



A. Stable modes 

Now the dispersion relation of the collective modes 
of the composite system is given by the zeros of Eqs. 
Jbl and (|3"6"|) . The propagator have poles at real value 



u > |k|. The dispersion relation for A-modes can be 
determined by finding the solution of the equation 



lo 2 a = k 2 + a(uj A ) + a'(uj A ) 
In case of G-modes we factorize A^ 1 as 



where 



and 



uj 2 g± = -(id 2 ±\Jn 2 + 4k 2 n 2 <5 2 ), 



Co 2 = a + a' + 13 + P' + 7 + k 2 , 
n = a + a' -P-P' + 7 + k 2 . 



(37) 



(38) 



(39) 



(40) 



For real uj > |k|, the square root of Eq. (|3^|) is always 
positive. Therefore, at most two stable modes come from 
G -modes. 

These solutions are found numerically by varying the 
various parameters introduced earlier. The results for the 
stable ^4-mode are shown in FigfT] for two values of the 
anisotropic parameter £ = {1, 10}, r\ — 0.2, 6j et = tt/2, 
Vjet = 0.7 and 9„ = with and without the jet. It is seen 
that the collective modes with the jet differs reasonably 
from that without the jet. Next we consider the depen- 
dence of the dispersion relation on the angle of propaga- 
tion of the jet with the wave vector. This is displayed in 
Figf2]for £ = 10, Vj e t = 0.7, 6 n = and 77 = 0.2. It is 
found that the dispersion relation is sensitive to 9j e t- In 
order to see the dependence of the collective modes on 
the jet velocity Vj e t, we plot the dispersion relation for 
A-mode for fixed 9 n and 9 jet in FigJ3] for two values of 



5 



1.25 




6 = 



e =ti/2 



FIG. 2: (Color online) Same as Fig[T]for £ = 10, v jet = 0.7, 
n = 0, rj = 0-2 and 6 jet = {0, tt/4, tt/2}. 
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FIG. 3: (Color online) Same as Fig[T]for v ]et = {0.7, 0.999}, 
£ = 10, 6 n = and 6 jet = tt/4. 



Vj^. The dispersion relation is modified again, but it is 
not that sensitive with Vj e t- The dispersion relation for 
the stable A-modes in the composite system is shown in 
the FigU]for 9 n = 0, and tt/2. We have fixed the param- 
eters as explained in the figure. It is observed that for a 
fixed value of 9 n the modes have significant dependence 
on the value of r\. When the wave vector is parallel to 
the direction of the anisotropy, i.e., k||n then 7 and n 2 
vanish identically. Therefore the dispersion relations for 
G+-modes are similar to A-modes for 9 n = 0. The de- 
pendence of the G+-modes on the anisotropic parameter 
£ is displayed in Figj5l For fixed 9 n we find significant 
difference in the collective modes with and without jet. 
However we do not find any significant dependence of the 
G+-modcs on the jet velocity (see in Fig|3]). From FigJTl 
it is clearly seen that the dispersion relation depends on 
the parameter r\ and the modes are strongly modified by 
this parameter. The collective modes also depend on the 
angle of propagation with respect to the anisotropy vec- 
tor, 9 n as depicted in Figj8] The numerical results for 
the stable G_-mode are shown in FigJS]for £ = {1, 10}, 
7] = 0.2, 9j et — tt/2, Vj et — 0.7 and 9 n = with and 
without jet. We find significant difference in the collcc- 



FIG. 4: (Color online) The dispersion relation for the stable 
A-mode for the system composed of anisotropic plasma and 
jet for £ = 10, v jet = 0.7, 8 jet = tt/2 and 77 = {0.02, 0.2}. 
The left(right) panel corresponds to 9„ = 0(n/2). 
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FIG. 5: (Color online) The dispersion relation for the stable 
G+-mode for an anisotropic plasma with jet for 9 n — tt/2. 
The left(right) panel corresponds to anisotropy parameter £ = 
1(10). 



tive modes with and without jet and with jet G^-mode 
quickly bends towards the lu = k line. The dispersion 
relation for G_ mode for 9j et = {0,7r/2} is displayed in 
FigHUl Unlike the A and G+-modes, it is seen that the 
dispersion relation for G_-mode is more sensitive to 9j e t- 
Like the A and G+-modes, the collective mode for G_- 
mode is also marginally sensitive to the jet velocity. In 
FigllTl it is observed that the dispersion relation for G_- 
mode strongly depends on the direction of the anisotropy 
axis. 
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FIG. 6: (Color online) Same as Fig[5]for v jei = {0.5, 0.999}, 
£ = 10, n = tt/3, Ojet = tt/2 and 77 = 0.2. 



FIG. 8: (Color online) Same as Figljfor 0„ = {0, tt/4, tt/2}, 
£ = 10, 7j jet = 0.7, 6> jet = 7r/2 and 77 = 0.2. 




k/m 



FIG. 7: (Color online) Same as FigE] for 77 = {0.02, 0.2}, 
£ = 10, 7j jet = 0.7, n = tt/2 and (9 jet = tt/2. 
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Unstable modes 



In the static limit, for small 9 n both the scale m 2 and 
m 2 _ are negative and so also a'. It indicates that the 
whole system is unstable with respect to magnetic insta- 
bility pjl [T3] - This can be identified as the so called 
filamentation or Weibel instability [HJ. The instability 
is driven by the energy transferred from the particles to 
the field. The growth rate of the filamentation instabil- 
ities is largest when the wave vector is along the beam 
line, i.e., k||n(0„ = 0) 0, MM 

in which case 7 and n 
vanish identically. Therefor the dispersion relations for 
the unstable modes are determined by the solutions of 
the following equations: 



cj 2 - k 2 - a{u) - a'{uj) = 
u? - p(u) - (3'{oj) = 



(41) 
(42) 



In the numerical simulation we find that unstable mode 
exits only for a — a' mode in the special case(k||n) con- 
sidered here. To solve the dispersion relation for the un- 
stable a — a' mode we first consider 9j et = and it is 
found that uj is purely imaginary i.e. uj = iT with T real 



FIG. 9: (Color online) The dispersion relation for the stable 
G_-mode for an anisotropy plasma and a jet for 9 n — tt/2. 
The left(right) panel corresponds to anisotropy parameter £ = 
1(10). 



valued. Results for different values of the anisotropy pa- 
rameter £ = {1, 10} are shown in FigfTJ] It is observed 
that instability first increases in comparison to the no 
jet case and then it becomes damped. To find the maxi- 
mum value of the momentum k max at which the unstable 
mode spectrum terminates, we take the limit T — > 0, to 
obtain(0jet = 0), 



v/?+(£-l)arctan(V?) 



"jet 



l D- 



(43) 



We shall now consider 0j e t dependence of the unstable 
modes. For nonzero 8j e t we find that collective modes are 
unstable propagating modes i.e. u> = a + iT. In Fig |13l 
we plot the imaginary part of the dispersion law of the 
unstable a — a'-mode for four different 9j et . It is clearly 
seen that T first increases with 9j et and then turns into 
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FIG. 12: (Color online) The growth rate F of the unstable a — 
a' mode, for two different anisotropy parameter £ = {1, 10} 
and 8j et = 0. 



FIG. 10: (Color online) The left(right) panel corresponds to 
the dispersion relation for the stable G_(G)-mode for the 
composite system for f9 jet = {0, 7r/2}, £ = 10, Vj e t = 0.7, 
9 n = 7r/3 and r\ = 0.2. 




FIG. 11: (Color online) Same as Fig|9]for £ = 10, v jet = 0.7, 
: tt/2, r) = 0.2 and different n = {0, tt/4, tt/2}. 



7 jet 



damping except 6j e t = tt/2 in which case the unstable 
mode always grows. It is interesting to see that when 
Ojet = t/2| the maximum value of T decreases and it 
never becomes damped. 

The growth rate dependence on the parameter r\ = 
{0.02, 0.2} corresponding to the value of the £ = 10 and 
9j et = is displayed in FigITU It is seen that instability 
grows with r\. 



V. SUMMARY 

We have calculated the dispersion relation of modes 
(both stable and unstable) due to the passage of rel- 
ativistic jets in an anisotropic plasma. It is seen that 
the dispersion relations of stable modes in an AQGP are 
modified due to the introduction of relativistic jet com- 
pared to the case when there is no jet. Dependences of 




FIG. 13: (Color online) Same as FigH2l with 9 jei 
{0, 7r/4, tt/3, tt/2}, ^ = 10 and t? = 0.2. 



possible modes on the direction of propagation of the jet 
and anisotropy direction have been demonstrated. We 
have also shown that the growth rate T increases with 
the passage of jet and this growth rate strongly depends 
on the anisotropy parameter £. This means that the in- 
troduction of jet in an AQGP leads to faster isotropiza- 




FIG. 14: (Color online) Same as Fig[TJ]with n = {0.02,0.2}, 
£ = 10 and 8 jet = 0. 
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tion for the special case considered here. It is also (3 — /3'-mode in the special case. However, it remains to 

interesting to note that the growth rate for 9j et = 7r/2 be seen how the unstable modes behave when a more 
never becomes damped. We also find no unstable for general case (arbitrary 8 n ) is considered. 
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